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ABSTRACT 

We present a Chandra study of the hot ISM m the giant eUiptical galaxy NGC4649. In common 
with other group-centred eUipticals, its temperature profile rises with radius in the outer parts of 
the galaxy, from ~0.7 keV at 2 kpc to ~0.9 keV by 20 kpc. However, within the central ^2 kpc 
the trend reverses and the temperature peaks at ~1.1 keV within the innermost 200 pc. Under the 
assumption of hydrostatic equilibrium, we demonstrate that the central temperature spike arises due 
to the gravitational influence of a quiescent central super-massive black hole. We constrain the black 
hole mass (Mbh) to (S.SSlo.gs) x lO^M© (90% confidence), in good agreement with stellar kinematics 
measurements. This is the first direct measurement of Mbh based on studies of hydrostatic X-ray 
emitting gas, which are sensitive to the most massive black holes, and is a crucial validation of 
both mass-determination techniques. This agreement clearly demonstrates the gas must be close to 
hydrostatic, even in the very centre of the galaxy, which is consistent with the lack of morphological 
disturbances in the X-ray image. NGC 4649 is now one of only a handful of galaxies for which Mbh has 
been measured by more than one method. At larger radii, we were able to decompose the gravitating 
mass profile into stellar and dark matter (DM) components. Unless one accounts for the DM, a 
standard Virial analysis of the stars dramatically over-estimates the stellar mass of the galaxy. We 
find the measured J-band stellar mass-to-light ratio, 1.37 ± O.IOMqLq^ , is in good agreement with 
simple stellar population model calculations for this object. 

Subject headings: Xrays: galaxies — galaxies: elliptical and lenticular, cD — galaxies: individual 
(NGC4649)— black hole physics 
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1. INTRODUCTION 



io boooming inoroopingly olopr that DupormaDaivo 



holes (SMBHs) in the centres of galaxies are in- 



t ima toly involved in the evolution of their hosts. Models 

of p;tJlaYy fnrmatinn gnggptit that pnnrgy injprtinn frnm 



active galactic nucleus (AGN) outbursts, in particular 
during hierarchical assembly, can have a signific ant im- 
pact on the structure of the evolv ing galaxy (e.g. ^ilk & 
Rees 1998; Di Matteo et al. 2005). In galaxy groups and 
clusters, cavities in the hot mtra-cluster medium (ICM) 
are a direct, observable manifestation of an episodic, 
AGN-driven feedb ack process acting on the ICM (e.g. 
Birzan et al. 2004). Such feed back has been implic ated 
in global heatmg of the gas (e.g. [ Donahue c £ al.| 2006|), re- 



distr i buting metals through the ICM (e.g. Mathews et al 

^004^ ) and 'Tnv n]<-pr\ a pnt;gihlp t;nliitinr i tn thp "rnni 



ing paradox" ( Mathews fc Brighenti||2003| , and references 
therein). 

Evidence for SMBHs has been found in increasingly 
large numbers of non-active galactic nuclei, leading to the 
picture of ubiquitous SMBHs in stellar sphero ids which 
has begun to emerge (Ferrarese & Ford [2005 , for a re- 
view). The black hole mass (Mbh) has been found to 
correlate with properties of the host, such as luminos- 



ity O f maoa of the bulge ( t4 ermondy fe Riohotone - 1095 
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199^ ; pVIcLure fc Dunlop] |2002| ; 



Marconi 



ic central light concentration ( Graham 



ct al. 2001), or the central velocity dispersion (a^,) of 



al 



Ferrarese & Mer- 



|2000a| ; 

ritt 2000| ; Tremaine et al.| ^002] ). These observational 
relationships are starting to bec ome a key constraint 
on galaxy-formation models (e.g. Robertson et al. 2006 ; 
Granato et al||2004| ). There remain however, concerns 
over the exact shape of the relations. In particular, dif- 



ferent au thors have found different slopes for the MBH-g ^, 
relation (Tremaine et al. 2002; Ferrarese & Ford 2005| ), 
and there are suggest ions it mav d eviate from a sim- 
ple power law relation ( Wyithc 2006 ). Possible disagree- 
ment in the masses of the largest SMBHs inferre d from 
the Mbh-c * and MBH-bulge luminosity relations (Lauer 



ct al 



^007 ) raise questions over which relation is more 
fundamental. Clearly, further progress requires a careful 
assessment of possible systematic biases in the black hole 
mass measurements, especially in the high-mass regime. 

A few independent techniques have evolved to measure 
the masses of SMBHs, each relying on its own set of sim- 
plifying assumptions. In galaxies with quiescent nuclei, 
the presence of the black hole is measured by its grav- 
itational infiuence on the kinematics of either the stars 
or, if present, a central ionized ga s disk (e.g. Gebhardt 



et al 



20031; pviacchetto et al. 
ound 



1997[ ). Since gas disks are 
stellar kinematical studies. 



not found in every galaxy, 
based on measuring the line-of-sight velocity dispersion 
profile are more generally applicable. Such studies are, 
however, complicated by a strong degeneracy between 
orbital structure and the gravitating mass profile, par- 
ticula rly in slowly-rotating galaxies ( Binney fc Mamon 
1982). Although measuring high-order moments of the 



2 



velocity dispersio n profile allows the degenera cy to be 
partially broken ( van der Marcl fc Franx , |1993[). th e so- 
lution is not, in general, unique uri et airp004|) . To 
obtain interesting constraints on the SMBH mass, there- 
fore, additional constraints such as a restricted range of 
allowed orb its or smoothlv-yarving phase-space are re - 
quired (e.g. Magorrian et al. 1998 ; Gebhardt et al. 2003 ) 



If th ese as sumptions are unrepresentative, however, Val- 



luri et al. pointed out that the systematic errors in the 
Mbh determination can exceed the statistical errors. 

Given the good agreement in the shape of the Mbh- 
fj, relation measured using masses obtained with dif- 
ferent methods (e.g. pebhardt et al. 2000b; Ferrarese 
et alj 2001), it seems unlikely that SMBH niasses Je^ 
term med from state-o f -the-art "3-integral" in ethods (van 
der 4larel et al.| |1998| ; [Gebhardt et~all |2003| ) are on av- 
erage in error. Nevertheless, a proper assessment of the 
systematics inherent in the mass determination, by com- 
paring masses for the same black hole measured using 
different techniques is of considerable importance. To 
date there are only a handful of galaxies for which re- 
cent ,«tp1lar Vinpmat.iral Mj^jj mpagin-pmpntg liavP hppn 



com pared with at least one reliable, independent tech- 
nique. Stellar and gas dynamical estimates of Mbh have 
been found to agree r easonablv in at least two svstems. 



^0061) 



NGC 3379 and Gen A ( [Shapiro et al.| p006[ ; [Marconi et al 



although there is evidence of non-circular motions 
in NGG3379, which can lead to a systematically mis- 
estim atcd Mbh. Simil ar agreement has been found for 
M81 ( Kormendy ^004 ), although the stellar kinematics 
measurement was based on "2-integral" models, which 
are potentially unreliable as they place constraints on 
the orbital structure which appear unrealistic, such as 
assuming the orbits are described entirely by only two 
inte grals o f motion (for more discussion see Ferrarese & 
Ford |200^ ) . Non-circular gas motions may explain strong 
discrepancies between Mbh measured by these diff erent 
tech niques in the g alaxies IG 1459 and NGC 4335 (Gap- 
^002r|Verdoes Klcijn et al.[ [2002[) 



pellad et al. 



this concern, the 



Despite 

bh-c* relation is usually constructed 



by combining samples of galaxies in which Mbh has 
been determined from stellar kinematics and samples for 
whic}i gas dynamics has been used (e.g. Tremainc et al. 
^002[ ), and there is little evidence that the sub-samples 
obey appreciably different relations. 

Other methods of determining Mbh are available in 
galaxies containing active galactic nuclei (AGNs), such 



as m easuring the motions of masing gas ( Miyoshi et al 
1995) or "reverber ation mapping" , usin g the variability 



of the AGN (e.g. Peterson et al. 2004). In two cases 



Mbh determined from reverberation mapping has been 



compared to stellar kinematics results ( Davies et al. 2006 
pnken et al. 2007 ), but since the uncertainties in the 
former techniqu e lead to masses a ccurate only to within 
a factor of ~3 ( pnken et al. 2004), these measurements 
do not provide good constramts on the validity of the 
stellar models. 

A new method to ident ifv SMBHs in giant e l liptica l 
galaxies was suggested by Brighenti fc Mathews (1999), 
who pointed out that the gravitational influence of a cen- 
tral black hole should have an impact on the distribution 
of the characteristically hot interstellar medium (ISM). 
The ISM should gradually flow inwards as it loses energy 
from the emission of thermal X-rays but, since the local 



TABLE 1 
Observation summary 



ObsID 


Start Date 


Exposure (ks) 


785 


2000 Apr 20 


18 


8182 


2007 Jan 31 


48 


8507 


2007 Feb 1 


15 


Total 




81 



Note. — Details of the Chandra observations used in 
the present analysis. For each dataset we quote the obser- 
vation ID number (ObslD), the start date and the exposure 
time, having removed periods of background "flaring" (Ex- 
posure) 

cooling time far exceeds the local free-fall timescale, the 
flow is highly subsonic and the gas remains close t o hy- 
drostatic equilibrium (Mathews & Brighenti ^003 ). In 



the powerful gravitational field close to the black hole, 
the inflowing gas should be sufficiently compressed to 
cause a strong temperature peak, which could be used 
as a black hole diagnostic. However, detections of cen- 



tral temperature peaks in galaxies are rare ( 


Humphrev 


et al. 


2004 ; Humphrey & Buote 


21 M Ml.,; 


Hump 


hrey et al. 



ified observationally. In part this reflects the difficulty 
in obtaining precise temperature constraints from X-ray 
spectroscopy on scales of a few hundred parsecs. AGN- 
induced disturbances are also common in the centres of 



early-type galaxies (Jones et al. 2002; Birzan et al. 2004) 
and these may disturb the gas enough to destroy many 
such hot cores. 

In this paper, we present a detailed Chandra study 
of NGG4649, a nearby, very relaxed and X-ray bright 
galaxy which is arguably the most promising candidate 
in which to sear ch for this effect . Base d on stellar kine- 
matical studies, Gebhardt et al. ( ^003 ) reported a black 
hole mass of {2.0toM) X IO^Mq (1 -a errors), placing it at 
the upper end of the Mbh-c* relation. A modest archival 
Chandra dataset has already revealed a remarkably re- 
laxed X-ray morphology (imp lving hydrostatic e quilib- 
rium is a good approximation: Buote & Tsai 1995) and a 
temperature profile which peak s in the centre ( Humphrey 



et al. 200e : RandaU et al. 2004). The object is faint in the 



radio ( Condon et al. 200^ ), indicating that it does not 
host a powerful AGN, and there is little evidence of mor- 
phological disturbances indicating AGN-ISM interaction 
(§ 2.1). On a large scale, excellent agreement between the 



mass determine d from globular clu ster kinematics and 
X-ray methods ( Bridges et al. 2006 ) also provides good 
support for hydrostatic equilibrium. We present deep 
new data which enable us to obtain interesting temper- 
ature constraints on scales as small as ^150 pc, which is 
crucial given the ~60 pc sphere of infiuence of the SMBH 
on the gas^. For the first time in any system we are able 
to place constraints on Mbh from hydrostatic X-ray gas. 

We adopted a nominal distance of 15.6 Mpc for 
NGG4649, based on the I-ban d SBF distance modulus 
estimate of Tonry et al. (2001), which was corrected by 
-0.16 to account tor recent revisions to the Gepheid zero- 

^ Defined as GMgu/c^, where Cg is the central gas sound-speed 
(in the absence of a black hole), ~ 490fcm from our data. 

This is slightly smaller than the sphere of influence on the stars, 



([Gebha 



ebhardt et al 



20o;) 



.^100 pc for the measured cr, = 385fcm, 
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Fig. 1. — Top left: J-band 2MASS image of NGC4649. Note the interloper spiral galaxy NGC4647 to the northwest. Top Right: point- 
source cleaned Chandra image of NGC4649, smoothed by a 5" Gaussian kernel, shown with arbitary, logarithmically-spaced contours. 
Note there is only a small amount of emission apparently associated with NGC 4647. Bottom left: "Unsharp-masked" Chandra image of 
the central part of NGC 4649 (note the different scale), created by subtracting an image smoothed with a 3" Gaussian kernel from one 
smoothed by a 1" kernel. The dark ring around the central, bright emission is an artefact of the smoothing process. Bottom right: "residual 
significance" image (see text) of the central part of the galaxy, indicating deviations from a model derived from elliptical isophote fitting to 
the data. The galaxy is very relaxed at all scales, and we do not confirm recent suggestions of morphological disturbances. Center: central 
22"x22" region of the residual significance image. 



point. At this distance, 1" corresponds to 76 pc. Unless 
otherwise stated, all error-bars represent 90% confidence 
limits. 

2. X-RAY DATA ANALYSIS 

The region of sky containing NGC 4649 has been im- 
aged by the Chandra ACIS instrument in ACIS-S con- 
figuration on three separate occasions, as listed in Ta- 
ble |ll We processed each dataset independently, using 
the (JIAO 3.4 and Heasoft 5.3.1 software suites, in con- 
junction with Chandra calibration database (Caldb) ver- 
sion 3.3.0.1. To ensure up-to-date calibration, all data 
were reprocessed from the "level 1" events files, fol- 



lowing the standard Chandra data-reduction threads^. 
We applied the standard correction to take account of 
the time-dependent gain-drift and charge transfer inef- 
ficiency, as implemented in the CI AO tools. To iden- 
tify periods of enhanced background ("flaring"), which 
seriously degrades the signal-to-noise (S/N) we accumu- 
lated background lightcurves for each dataset from low 
surface-brightness regions of the active chips, excluding 
obvious point-sources. Periods of flaring were identified 
by eye and excised. The final exposure times are listed 

^ http:/ /cxc. harvard.edu/ciao/threads/index. html 
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in Table 0. 



2.1. X-ray image 

For each dataset we generated a full resolution im- 
age in the 0.3-7.0 keV energy-band and a correspond- 
ing exposure map computed at an energy of 1.7 keV. 
Point sources were detected in each image with the CIA 
wavdetect task, which was set to search for structure at 
scales of 1, 2, 4, 8 and 16 pixels, and supplied with the 
exposure-maps to minimize spurious detections at the 
image boundaries. The detection threshold was set to 
10~^, corresponding to < 1 spurious source detections 
per chip. To improve the signal-to-noise (S/N) of the 
data, we combined the resulting images and exposure 
maps. Given errors in the absolute astrometry of Chan- 
dra, we corrected the aspect onto a common frame of 
reference before addition. To achieve this registration, 
we matched common point-sources in each image and 
shifted the data to ensure the maximum number of such 
sources overlapped. Since there can be variability in a 
significant fraction of the LMXBs in any galaxy (or, in- 
deed, in background AGN) we did not require all point 
sources to match. Having generated a merged image and 
exposure-map, to improve signal-to-noise (S/N), we re- 
detected the point sources. All detected sources were 
confirmed by visual inspection, and, for each, appropri- 
ate elliptical regions containing approximately 99% of its 
photons were generated. 

We examined the image for evidence of morphological 
disturbances which may suggest deviations from hydro- 
static equilibrium. First we processed each image to re- 
move the point-sources. To do this, we replaced all pho- 
tons within the source detection ellipse of each source 
with a rtificial dat a , usin g the algorithm described in de- 
tail in Fang ct al. (2008). In short, since the wavdetect 
task can be used to generate a smoothed, "normalized" 
source-free image, we replaced the data within each 
source detection ellipse from our real image with data 
from the corresponding pixels on the source-free image, 
after adding Poisson noise and taking exposure-map vari- 
ations into account. We found this procedure generally 
worked extremely well at removing sources. The data 
were then fiat-fielded with the exposure map and exam- 
ined by eye. We show in Fig |l the image of NGC4649, 
having smoothed it slightly with a 5" (0.4 kpc) Gaussian 
kernel. The image is remarkably round and shows no 
evidence of disturbances or structure. It is interesting 
to note that the X-ray image is rounder than the optical 
light (which is, nonetheless, quite round, with a major/ 
minor axis ratio of < 1.2). This is unsurprising since gas 
in hydrostatic equili brium traces the potentia l (the "X- 
ray shape theorem" : Buote fc Canizare^ 1994 ) , which is 
rounder than the mass distribution giving rise to it. 

To seek evidence of more subtle structure, we exper- 
imented with "unsharp masking" the image, i.e. differ- 
encing the image smoothed at two different scales. We 
did not find convincing evidence of any morphological 
disturbances for any of the scales (using smoothing ker- 
nels from 3-20") we tried. We show in Fig |l] a typical 
"unsharp-masked" image, illustrating the lack of sub- 
structure (note that the dark ring around the central 
bright region i s an artefact of the smoothing process). 
^hurkin et al. (200S) recently argued for the presence of 
AGN-induced cavities in the centre of NGC4649, but. 



using deeper data than these authors, we found little ev- 
idence for such features. Performing the "unsharp mask- 
ing" they describe, we found little evidence for the sup- 
posed cavities, save features which are practically indis- 
tinguishable from erratic structure arising as an artefact 
of the image processing. Nonetheless, to assess whether 
any such depressions do exist, we first fitted elliptical 
isophotes to the central part of the X-ray image with the 
PC-IRAF task ellipse. We then extracted the total 
counts both in the image and the model within two circu- 
lar apertures of 4" di ameter placed appro ximatelv at the 
positions reported by ^hurkin ct al. ( 2008| ) for the depres- 
sions. One can then compute =(data-model)^/model 
which can be used to assess the significance of the fea- 
tures. A standard would assume that we had per- 
formed a random trial, however, which is clearly incor- 
rect here. Instead, we performed Monte-Carlo simula- 
tions to assess the true distribution of . We performed 
1000 such simulations, in each of which an artificial im- 
age was generated by adding Poisson noise to the model 
image. The model was then subtracted and the resulting 
image smoothed (with a 2" Gaussian) to bring out low 
counts structure at the scale of interest which lay close to 
the apparent north and south radio lobes (by which we 
mean within 25"x25" boxes centred on each lobe). We 
repositioned the circular apertures at the corresponding 
minima and computed from the (unsmoothed) artifi- 
cial data. We found that 5.3% of the simulations gave 
larger than we obtained with the real data, implying 
that the putative cavities are significant only at 
We therefore consider the evidence for AGN-induced cav- 
ities in the hot gas to be weak. 

To investigate the presence of other possible structures, 
we show in Fig |^ an indication of deviations from our 
model fit by plotting (data-model)^ /model, which cor- 
responds to the x^ residual in each pixel. To bring out 
structure, we smoothed this image with a Gaussian ker- 
nel of width 0.5" (1 pixel). As is clear from this "residuals 
significance" image, there is no evidence of any coherent 
residual features (with the possible exception of one or 
two point-sources that may have been imperfectly re- 
moved). The X-ray emission is therefore very smooth 
and symmetric. ^ ^ 

Randall et al.' ( 200(^ , p004D adaptively smoothed Chan- 
dra and XMM images ol NGC 4649 (using shallower 
Chandra data than the present work) and identified ra- 
dial features they dubbed as "fingers of emission" . How- 
ever, although adaptive smoothing can give a useful over- 
all impression of the data, it tends to introduce artificial 
structure into images. Our data are deeper and we do not 
find convincing evidence of such features. Furthermore, 
visual comparison of the authors' Chandra and AMM im- 
ages suggests some differences; in particular the radially- 
averaged, azimuthal surface-brightness plots do not obvi- 
ously resemble each other, with the AMM data appearing 
(by eye) to show a larger amplitude variation with po- 
sition angle than the Chandra data. Since XMM is less 
able to resolve out LMXBs, we speculate that unresolved 
point-sources may be largely responsible for the signifi- 
cant structure reported for the XMM data. Even taking 
the supposed XMM structures at face value, they would 
imply azimuthal variations in gas density of only'^5%. 
Small variations of this magnitude do not constitute sig- 
nificant deviations from hydrostatic equilibrium and are 
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smaller, in general, than our r neasurement er r ors. B ased 
on the shallower observation, Randall et al.| ( ^004 ) also 
report a ~ 2.5-a detection of a central ~5" structure, but 
inspection of the "residual significance" image shown in 
Fig[^ indicates there are no significant residual structures 
in the central region. 

2.2. Spectral analysis 

We extracted spectra in a series of concentric, con- 
tiguous annuli, placed at the X-ray centroid. We deter- 
mined the centroid iteratively by placing a 0.5' radius 
aperture at the nominal galaxy position (obtained from 
NED) and computing the X-ray centroid within it. The 
aperture was moved to the newly-computed centroid, and 
the procedure repeated until the computed position con- 
verged. The widths of the annuli were chosen so as to 
contain approximately the same number of background- 
subtracted photons, and ensure there were sufficient pho- 
tons to perform useful spectral-fitting (in this case we had 
> 2000 background-subtracted photons per spectrum) . 
We placed a lower limit of 2.5" on the annulus width, to 
ensure that the instrumental spatial resolution does not 
lead to strong mixing between the spectra in adjacent an- 
nuli. The data in the vicinity of any detected point source 
were excluded, as were the data from the vicinity of chip 
gaps, where the instrumental response may be uncertain. 
To prevent possible low-level contamination from any hot 
gas in the central part of the interloper galaxy NGC 4647, 
we additionally excluded data in its vicinity from the 
9.3-13.1 kpc annulus, which intersects its centre. We ex- 
tracted products from all active chips, excluding the S4 
(which suffers from considerable "streaking" noise). Ap- 
propriate count- weighted spectral response matrices were 
generated for each annulus using the standard CI AO 
tasks mkwarf and mkacisrmf . For each spectrum, we 
estimated the ba c kgrou nd using the method outlined in 



Humphrey et alj (2006). We extracted identical prod- 
ucts individually for each dataset, taking into account 
the astrometric offsets between them (determined from 
the image registration discussed above). The source and 
background spectra were added, and the response ma- 
trices averaged using the standard Heasoft tasks addrmf 
and addarf , with weights based on the number of pho- 
tons in each spectrum. We found there was very little 
emission outside ~5' (23 kpc), indicating truncation of 
the X-ray halo outside this radius as NGC 4649 falls into 
the Virgo cluster. 

Spectral- fitting was carried out in the energy-band 0.5- 
7.0 keV. The spectra were rebinned to ensure a S/N 
ratio of at least 3 and a minimum of 20 photons per 
bin (to validate fitting). We fitted data from all 
annuli simultaneously using Xspec. To model the hot 
gas we adopted a vapec component, plus a thermal 
bremsstrahlung component to account for undetected 
point-sources (this model gives a good fit to the compos- 
ite spectrum of th e detected sources in nearby galaxies: 



[rwin et al. 


^0031). In common with our orevious 


anal 


ysis (|Hump 


tirey & Buotc 


L'iMM,.,: 


Humphrey et al. 


2004 



we adopted a slightly modified form of the Xspec vapec 
code to enable us to tie the ratio of each elemental abun- 
dance with respect to Fe between each annulus, while 
the Fe abundance (Zpc) was allowed to fit freely. To im- 
prove S/N Zpc was tied between adjacent annuli where 
constraints in individual annuli were poor. We allowed 



the global ratios of O, Ne, Mg, Si, S, and Ni with re- 
spect to Fe to fit fr eely, and fixed the r emaining ratios 
at the Solar value ( Asplund et al. 2004). The absorb- 
i ng column density (Nh ) was fixed at the Galactic value 
( Dickey fc Lockniai^ 199C| ). To account for projection ef- 
fects, we used the projct model implemented in Xspec. 
Since the X-ray emission appears truncated at large ra- 
dius, we did not include a component to account for pro- 
jection from outside our outermost annulus. However, as 
NGC 4649 is within the Virgo cluster, we took account 
of possible interloper cluster emission by including an 
addi tional hot gas component, w ith kT fixed at 2.5 keV 
(e.g. [Gastaldello fc Molendi|[2002| ) , which was assumed to 
have a constant surface brightness over the field of view. 
We obtained a reasonable fit (x^/dof=1178/1088) to the 
spectra with this model, which is of co mparable quality 
to ou r previous analysis of 7 systems ( prlumphrey et aL 
2006). The best- fitting abundances were in agreement 



with those reported for NGC 4649 in that paper. Error- 
bars we re computed via t h e Mo nte-Carlo technique out- 



lined in Humphrey et al. (E006), and we carried out 50 



error simulations. The resulting temperature and density 
profiles are discussed below. 



3. MASS MODELLING 



3.1. 



In 



Entropy- Temperature method 
our previous work ([Humphrey et al. 2006 



Gastaldello et al. 2007; ^appacostact al. 2006) we consid- 
ered a variety of difterent approaches to solve the equa- 
tion of hydrostatic equilibrium and obtain the mass pro- 
file. In our past analysis of NGC 4649, we adopted a pa- 
rameterized model for the total gravitating mass profile 
and an ad hoc model for the temperature profile, solving 
for the total gas density. Although this method allowed 
us to obtain interesting constraints, it relies to some ex- 
tent on the adoption of a temperature profile with rea- 
sonable asymptotic behaviour. We found a wide range of 
different temperature profiles for the galaxy-scale objects 
in our sample, unlike the approximately "universal" tem- 
perature profile s of galaxv clusters and, to some extent 
galaxy groups (Vikhlinin et al. ^005 ; Gastaldello et al 
2007| ; [Rasmussen fc Ponman] |2007D 



"Xs an alternative to these methods, in the present work 
we chose to exploit the fact that the entropy profile is, 
in general, far smoother than the temperature or density 
profiles and can typically be fitted by simple parameter- 
izcd models which a rc well-behaved asymptotically (e.g. 
Donahue et al, ^006 ). In what follows we define the quan- 
tity S — p^^^^kT/l^mH), where p is the gas density, T 
is the temperature, k is Boltzmann's constant, mn is the 
mass of a hydrogen atom and fi is the mean molecular 
weight factor, taken to be 0.62. This is related to the 
actual entropy by a logarithm and multiplicative con- 
stant and is proportional to the usual "entropy" proxy 

used in the literature, s = Ue ^^^kT (e.g. Sanderson et al. 
2003 ). An advantage of considering the entropy profile 



is that stability against convection requires dS/dr to be 
positive at all radii, if the hydrostatic equilibrium equa- 
tion is valid, placing a further constraint on meaningful 
models and improving our mass constraints. 
One can write the equation of hydrostatic equilibrium 
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in terms of pressure, P, and S as 

dx^_2GMTOT(r)g-3/5 
dr 5 r^ 

where x—P^^^, G is the universal gravitating constant, 
and A/tot (r) is the total mass enclosed within a given 
radius r. For a given distribution of S and Mtot, this 
equation can be directly solved for x, provided the mass 
contribution of the gas can be neglected at all radii. Al- 
though this is true in the present work, since the gas 
contributes < 1% of the mass at all radii, in practice 
we allowed for a nonzero gas mass by differentiating and 
rearranging the equation to give: 



dr\ drj^ 553/5 ^ 5^ dr 



(2) 



where the mass enclosed, M, comprises only the dark 
matter, stellar mass and the black hole mass. The gas 
mass is represented by the second term on the left hand 
side. For any parameterized distributions of S and M, 
this equation can be solved uniquely for x using a Runge- 
Kutta method. This means that for any such param- 
eterized models of mass and entropy, the temperature 
and density profiles can be uniquely determined. We 
began all integrations at 10 pc, assuming the gas con- 
tributes negligibly to the total mass within this radius. 
The starting value for x (i.e. central pressure) at 10 pc 
was determined as a free parameter of our fit. 

For each annulus in which we extracted data, we explic- 
itly integrated appropriately weighted functions of pres- 
sure and temperature to compute the predicted contri- 
bution to the emission measure (and hence mean den- 
sity) from gas in the corresponding shell (recalling that 
the data have been deprojected). We also computed an 
emission-weighted temperature (ignoring the tempera- 
ture and abundance dependence of the emissivity, since 
neither term varies much over the field of view). We com- 
pared these averaged quantities to o ur measured tem- 
peratu re and density data-points. See iGastaldello ct al 



(2007, Appendix B) for a detailed discussion of incorpo- 
rating the plasma emissivity into the gas modelling. 

3.2. The stellar mass component 

In order to place constraints on the black hole mass 
in the centre of a galaxy, it is imperative that we have 
an accurate picture of the stellar mass in that region. 
To do this, we assumed that the stellar mass follows 
the three dimensional distribution of the optical light, 
provided we adopt data observed in as red a filter as 



possi ble. From our Lick index analysis (Humphrey et al 



^0061) 



we were able to estimate a mean age for the stellar 
population, as well as both a central stellar abundance 
and an emission- weighted abundance for the gala xy as a 



who l e. Ad opting the stellar population models of Maras- 
( 2005| ), the J-band M/Lj ratio should only change 



by ~3% depending on which of the two abundances is 
assumed. This supports the assumption that mass fol- 
lows light. To ensure maximum accuracy at the scale 
of our innermost Chandra data-bins, we first considered 
the HST WFPC2 data with the F814W (approximately 
I-band) filter (fo r more details on the HST analysis see 
Humphrey 2007 ). We considered only data on the PC 
chip (to prevent the particular geometry of the WFPC2 



plane complicating our analysis), limiting our profile to 
a radius of ~11". To extend the coverage to larger radii 
(^4'), we adopted the publicly-available 2MASS J-band 
images, masking out all data, with in the D25 ellipse (as 
listed in de Vaucouleurs et _al. 1991) of the interloper spi- 
ral galaxy NGC 4647 (Fig as well as from the vicinity 
of any bright point-sources. Although we were not able 
exactly to match the filters used in both cases, when the 
I-band HST profile is suitably rescaled, the profiles agree 
well outside ~2", implying that there is not a strong I- J 
colour gradient over the HST field of view, and hence we 
can use the HST data as a proxy for a high-resolution 
J-band image. 

By visual inspection, the axis ratio (major/ minor) 
of the optical isophotes is <1.2 at all radii of interest 
(Fig indicating that the galaxy is fairly round; in 
what follows we adopt a spherical approximation. First 
we investigated whether simple analytical models can 
adequately match the observed surface brightness pro- 
file. We found that the de Vaucouleurs model signifi- 
cantly over-estimated the optical light in the centre of 
the galaxy. A Sersic profile with n~2.6 was able to char- 
acterize the overall shape of the profile but, especially in 
the inner regions, could not capture its fine detail (Fig|^). 

Rather than using these simple presciptions, we instead 
adopted a numerical deprojection technique, which is jus- 
tified given the high quality of the data. We divided the 
emission into a series of concentric shells with outer and 
inner radii corresponding to the annuli in which our sur- 
face brightness profile (Fig ^ was computed. In each 
shell, we assumed the luminosity density of the stars 
was constant, enabling us to relate the observed surface 
brightness in each annulus to an appropriately weighted 
sum of emissio n from the shell i t intersects and those 
outside it (e.g. Kriss et al. 1983| ). We assumed there 



was no emission outside the outermost shell, allowing 
us to compute the luminosity density (in Lq kpc^^) in 
the last shell directly from the measured surface bright- 
ness in our outer annulus. We then stepped inwards, 
iteratively solving for the luminosity in each shell. This 
technique is effectively the same as the "onion p eeling" 
algor ithm widely used to deproject spectra (e.g. Buote 
2000). We assumed a J-band solar magnitude of 3.71, 
or consis t encv with the stellar population models of 
Maraston ( 2005 ). We have confirmed that the integrated 
ight profile is in good agre ement with the V -band depro- 
jected light profile used by Gcbhardt ct al. f 2003 ) (when 
scaled by a factor 3.34 which, using Maraston 's mod- 
els, accounts for the different photometric bands used 
given the metallicity and age of the stars in the centre 
of NGC 4649). The profiles differed by less than 20% 
over all scales at which they overlapped, and typically 
agreed within a few percent (especially in the crucial in- 
ner kpc). Such agreement provides further support for 
our assumption that the stellar M/Lj ratio is approxi- 
mately constant with radius. 

4. RESULTS 

From our spectral-fitting results, we have determined 
mean temperature and density data-points within each 
annulus. We converted these data to a mean (emission- 
weighted) entropy profile, as shown in Fig^. We chose to 
fit this profile with a simple empirical model, for which 
we found a broken-powerlaw plus a constant model to be 
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TABLE 2 

The best-fit results 



Parameter 



Value 



18.3/20 

1.37±0.10MoLq1 
2.00to 37 keV cm2 
SStfl keV cm2 
1.5 ± 0.18 
2A4tlZ kpc 
1.07±0.05 



Fig. 2. — Le/t.' Surface brightness profile of the galaxy, in optical light, shown with the best-fitting Sersic model. The Sersic model 
describes the overall shape, but not the detail of the luminosity profile. To aid the fitting, a systematic error of 5% has been added in 
quadrature with the statistical errors on each data-point. The J-band 2MASS profile is shown in black, and the renormalized I-band HST 
data in red. Right: The non-parametrically deprojccted J-band luminosity density profile. 

ratio, central black hole mass (Mbh), and total mass 
and concentration of the NFW dark halo to fit freely. 
We find that the global parameters of the dark mat- 
ter (NFW) halo are highly degenerate with each other 
without the applicati on of additional constraints (e.g. 
Hmnphrcy ct al. 2006 ), but their exact values have little 
bearing on the shape of the dark matter profile in the 
region of interest, so we will discuss them in detail in a 
future paper. 

Because the measured entropy and temperature in each 
annulus were, by definition, averaged quantities, in con- 
trast to the continuous, parameterized functions being 
fitted, care was taken over comparing the model with 
the data. The model density profile (squared) was ex- 
plicitly integrated over each bin to compute the emission 
measure, while the temperature model was similarly inte- 
grated with Pg^g (emission) weighting. These quantities 
correspond to the observables determined directly from 
the spectral fitting, and so could be used to compute 
the emission-weighted entropy (and temperature) for di- 
rect comparison with the data. The best-fitting entropy 
model is shown in Fig |^, and exhibits asymptotic (large 
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stellar M/Lj 
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SI 

Pi 

^brk 
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Note. — Selected best-fitting 
rameters (see text). 



pa- 



sufficient, i.e.: 



So + si/(r) 



(3) 



where: 



\ Sbrk 



{r < Sbrk) 
{Sbrk < r) 



(4) 



As explained in § pi for a pair of parameterized mod- 
els for entropy and the mass profile, one can uniquely 
determine the temperature (and density) profile. It is 
therefore possible to determine the best-fitting parame- 
ters of the mass model by fitting the entropy and tem- 
perature data simultaneously, while varying the parame- 
ters which describe the entropy and mass profiles. Since 
the fractional error-bars on the density are typically far 
larger than those on the temperature (due to partial de- 
generacies between the continuum and the unresolved 
point-source component, and between the gas density 
and abundance, for the range of temperatures considered 
here), we found this was more constraining than fitting, 
for example, entropy and pressure at the same time. 

To model the mass profile, in addition to the mass of 



scribed in Si 


3.2 




ing |Navarro 


et 


al. 



|I997| ), and a central supermassive 
black hole. We allowed the stellar mass-to-light (M/Lj) 



radius) behaviour consistent with the expected s oc 
from structure formation r nodels incorporating p urely 
gravitational processes (e.g. Tozzi fc Norman 2001). We 
also show the temperature and gas density profiles of the 
object, along with the best-fitting models in Figs ^ and 
Fig|. 

The radial dependence of the best-fitting mass model 
is shown in Fig ^, along with the contribution of each 
of the separate mass components. We find that the stel- 
lar mass component dominates the mass profile within 
~6 kpc, with a significant contribution of dark matter 
required to explain our observations at large radii. This 
is largely consistent with our previous analysis, despite 
the much simpler stellar mass profile we adopted in that 
work. Within our innermost annulus, however, we found 
that the rapidly- falling stellar mass was unable to explain 
the strong peak in the temperature profile. We show in 
Fig ^ the best-fit temperature model with the black hole 
mass (Mbh) fixed at zero, which clearly does not repro- 
duce the sharp central spike and is only a marginally ac- 
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Fig. 3. — The entropy profile of the data, shown along with the 
best-fitting model. Although the profile fiattcns at small radii the 
entropy gradient is always positive, as expected for gas in hydro- 
static equilibrium. 
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Fig. 4. — Measured temperature profile of the hot gas in 
NGC4649 (error diamonds). Also shown are the best-fitting mod- 
elled profile from the entropy-temperature analysis (solid line; 
black) and the best-fitting temperature profile if a central black 
hole is omitted (dashed line; red). In addition, the dotted line 
(hliuf ) indicates the best-fitting arbitrary parameterized model (See 

ceptable fit (global x^/dof=33.7/21). Allowing the black 
hole mass to fit freely enabled us to obtain an excellent 
fit to the data (x^/dof=18.3/20), an improvement sig- 
nificant at 99.94% significance, based on an f-test. The 
best-fitting values of key fit parameters are shown in Ta- 
ble |. 

4.1. Traditional mass analysis 

In our analysis, we adopted a new approach to X-ray 
mass analysis, which we expect to give a more precise de- 
termination of Mbh than would be expected from other, 
more traditional approaches. However, it is also of in- 
terest to determine whether interesting results can be 
obtained using the older methods, which have the ad- 
vantage of less computational expense. In order to in- 
vestigate this, we additionally analysed the data using 
a more orthodox mass-analysis approach (here dubbed 
the "traditional method"). In this case, we adopted ad 



Fig. 5. — Measured gas density profile of the hot gas in NGC 4649 
(error diamonds). Also shown are the best-fitting modelled profile 
from the entropy-temperature analysis (solid line; black) and the 
begt-fitt.ing arbitrary parameterized model (dotted line; blue. See 
§ 4.1). 
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Fig. 6. — Mass profile of NGC 4649 determined from our fitting. 
The individual data-points were.dpri^^pd from thr. traditional non- 
parametric method discussed in Gastaldello et al. (E007). We also 
show our best-fitting mass-profile determmed trom our entropy- 
temperature modelling method (see text). We show separately the 
total mass enclose (solid line; black), the stellar component (dashed 
line; red), the dark matter (dotted line; blue), the gas (dash-dot- 
dot-dot line; magenta) and the black hole mass (dash-dot line; 
orange). We stress that the models do not represent a fit of the 
data-points, which were derived independently. 

hoc parameterized models for both the temperature and 
density profiles. These models were chosen to capture 
as much detail of the profiles as possible, but may not 
be a perfect fit to the data. One can simply compute 
the mass enclosed within any radius by differentiating 
these models and inserting the derivatives into an appro- 
priately re-arranged form of the hydrostatic equilibrium 
equation. We used this approach to generate mass "data- 
points" corresponding to each of our annuli. We discuss 
this met hod, and the mod e ls we used to fit our data, in 
detail in [Gastaldello ct al.| ( |2007| ). 

For our present analysis, we fitted the density pro- 
file with a standard beta-model functional form, with an 
added central cusp. To fit the complicated temperature 
profile, we adopted a model comprising two smoothly- 



9 



TABLE 3 

Systematic error budget 



Tost 


Mbh 


M/Lj 




(IO^Mq) 




Best-fit 


3.35 


1.37 


AStatistical 


+0.67 
-0.95 


±0.10 


AMethod 


+0.75 


-0.13 


ADistance 


+0.01 
-0.27 


+0.06 
-0.12 


APlasma code 


+0.11 


-0.07 


ABandwidth 


-0.61 


+0.06 


AXRB 


0.005 


0.001 


ANh 


0.001 


0.001 


ABackground 


-0.10 


+0.22 


ACentroiding 


0.001 


0.02 


AAsymmetry 


+0.16 


-0.002 



Note. — Summary of the error-budget for the mass 
analysis. We Ust the best-fitting values of Mbh and 
mass-to-light (M/Lj) and the (90%) statistical errors 
(AStatistical). In addition we show estimates of the mag- 
nitude of the systematic uncertainties due to various fac- 
tors. These should not be added in quadrature with the 
statistical errors. We show the impact of adopting differ- 
ent fitting methods (AMethod; see text), the estimated 
errors due to distance errors (ADistance), uncertainties in 
the plasma code (APlasma code), the choice of bandwidth 
(ABandwidth), modelling of unresolved X-ray binaries 
(AXRB), the adopted Nh(ANh), treatment of the back- 
ground (ABackground), centroiding errors (ACentroiding), 
and possible asymmetry in the data (AAsymmetry). 

joined powerlaws. We show the fits to each profile in 
Figs ^ and |[ and the mass data-points are overlaid on 
our best-fitting mass model in Fig pi Clearly, there is ex- 
cellent agreement between these two profiles, indicating 
the robustness of our mass analysis results. There are, 
nonetheless, slight differences if we fit the mass data- 
points directly with our parameterized mass model (we 
omitted the gas mass, which is negligibly small over the 
range of interest), in large part because our adopted tem- 
perature profile does not fully capture the curvature of 
the data around 1-2 kpc. Nonetheless, from these data 
we found Mbh=(4.1 ± 2.2) x lO^Mg, which, although 
slightly higher than our best estimate, is statistically 
consistent with it. It is also within 1.6-c r of the stellar 

(|2003|) . 



kinematical estimate of [Gebhardt et al.| ( |2003|) . Simi- 
larly, our measured M/Lj ratio, 1.24 ± 0.22MqLq^ was 
slightly discrepant, but statistically consistent, with our 
preferred value. 

5. SYSTEMATIC ERROR BUDGET 

In this section we address the sensitivity of our results 
to various data-analysis choices which were made. We fo- 
cused on those systematic effects we estimate to be most 
likely to affect our results. An estimated upper limit on 
the sensitivity of our results to these choices is shown 
in Table |3[ These numbers reflect the sensitivity in the 
best-fit parameter to each potential source of systematic 
error, and we stress they should not be added in quadra- 
ture with the statistical errors. We outline below how 
each test was performed. Those readers uninterested in 
the technical details of our analysis may wish to proceed 
directly to § ^. 

First of all, in Table ||, we list (AMethod) the impact 
of u sing the traditional mass analysis method described 



details of the temperature and density profiles. Nonethe- 
less, the results are in good agreement from both tech- 
niques. Next we investigated the impact of changing the 
distance by ±1.8 Mpc, corresponding approximately to 
the 90% sta tistical error the I -band SBF distance mea- 
surement of Tonry et alj ( |2001| ) . This ran ge encompasses 
the ad opted distance (16.8 Mpc) used by Gebhardt et al. 
(2003). This had only a minimal impact on the measured 



Mbh, although more effect on the M/Lj ratio 

We next considered whether our spectral-fitting choices 
could have led to a systematic mis-estimation of the 
density and temperature, and hence Mbh- There are 
some uncertainties in the modelling of individual emis- 
sion lines, in particular for the important Fe blend, so we 
therefore experimented with replacing the APEC plasma 
model with the MEKAL code, which models the atomic 
physics differently. This had only a slight impact on 
the best-fit results (APlasma code). To assess whether 
the adopted bandwidth can affect our results, we exper- 
imented with fitting the data in the ranges 0.7-7.0, 0.5- 
2.0 keV and 0.4-7.0 keV, in addition to the default 0.5- 
7.0 keV. We estimated the impact on our results of im- 
properly modelling the unresolved X-ray binary spectral 
component by varying the temperature of the parame- 
terized (bremsstrahlung) model by ±25%, or replacing 
it with a powerlaw with r=1.5. Finally, to determine 
whether systematic errors in the adopted absorbing col- 
umn could lead to a large error on Mbh. we perturbed 



Nh b v ±25% from its nominal value ( Dickey • 
19901) . 



Lockman 



4.1 



as opposed to the default method we adopted. 
The traditional method is limited in that the adopted pa- 
rameterized models have difficulty in capturing the fine 



A potentially serious cause for concern when perform- 
ing spectral analysis of diffuse emission is the treatment 
of the background. We therefore investigated whether 
our choice for the background made an impact on our 
results by ad opting, instead of our preferred treatment 
given in § |2.2| , spectra extracted from appropriately pro- 
jected background template events files distributed with 
the Caldb. These template files are constructed from 
processed, stacked observations of nominally blank-sky 
fields. Although they do not directly take into account 
the large-scale variation in the cosmic X-ray background 
(nor variability in the non X-ray background compo- 
nent), they are widely used to estimate the background 
in Chandra analysis. The templates were renormalized 
to ensure the 9.5-12 keV count-rates of the source and 
the background agreed; this is formally incorrect since 
it also renormalizes the cosmic X-ray background in the 
templates as well. Nonetheless, we found that the choice 
of background did not have a significant impact on Mbh, 
although the systematic effect on M/Lj is larger than the 
statistical errors. This is easily understood since Mbh is 
determined by data from the bright, central regions of 
the galaxy, where the background is negligible. 

Another possible cause of systematic error in our anal- 
ysis is improper centroiding; if our spectral extraction re- 
gions were offset from the dynamical centre of the galaxy, 
this will cause the derived temperature and density pro- 
files to be systematically incorrect, leading to an error 
in the computation of the central mass. In practice, 
we would expect a mis-alignment to result in a system- 
atically under-estimated central temperature, reducing 
the significance of the SMBH detection. We checked the 
agreement of the Chandra centroid and the optical cen- 
tre (identified by eye from the HST images); based on 
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a bootstrap alignment of the Chandra and HST images, 
which r equired a fraction of the poi nt sources in each to 
match ( Humphrey fc Buote 2GG6b| ), we estimated they 
are consistent within ^0.6", which is comparable to the 
on-axis PSF of Chandra. We found that moving the X- 
ray centre by this amount made little difference to our 
overall result (ACentroiding), which is unsurprising since 
the shift is very much smaller than the radius of the cen- 
tral extraction region (2.5"). 

Finally, we considered whether low-level asymmetries 
in the X-ray image might constitute a strong deviation 
from our assumption of spherical symmetry. To investi- 
gate the possible impact of such features, we artificially 
restricted the field of view, re-extracting all spectra from 
a set of suitably-oriented concentric semi-annuli. We 
found the fit results agreed very well with our best-fit, in- 
dicating that any such asymmetries, if present, only have 
a minor impact. It is worth noting that, if the system 
globally deviated strongly from hydrostatic equilibrium, 
it would be highly unlikely that we would observe such 
symmetry. 

6. DISCUSSION 

Using the assumptions of hydrostatic equilibrium, 
spherical symmetry, a single-phase ISM and that the stel- 
lar mass follows the optical light, we have obtained an 
estimate for the mass of the supermassive black hole at 
the centre of NGC4649'^. This is the first time that the 
mass of a supermassive black hole has been determined 
from the distribution of hydrostatic, X-ray emitting gas. 
In principle, this method can also be applied to other gi- 
ant elliptical galaxies with relaxed X-ray morphologies, 
allowing the most massive black holes to be studied. 

We placed a lower limit on the black hole mass of 
1.6 X 1O^M0, at 99% confidence. Our measured mass 
agrea o (within " ■' 2 cr) with that determined from otel — 
lar Mnematical atudiea, (2.0tQ g) x 10^M(j (1 er error), 
based on a different set of assumptions. This makes 
NGC4649 one of only a handful of galaxies for which 
Mbh has been measured via multiple complementary 
methods. Although the overall agreement in the shape 
of the Mbh-c* relation when measured using different 
techniques (e.g. Gebhardt et al. |2000b| ) provides good 
evidence that the masses from stellar kinematics are not 
on average in error, the systematic uncertainties on any 
individual measurement of Mbh are more difficult to as- 
sess. In the few systems where the mass can be inferred 
from both gas and stellar dynamics, the results have been 
mixed, possibly on account of non-circular gas motions 
(§ 0). Therefore, the good agreement we found between 
the measurements of Mbh from stellar dynamics and hy- 
drostatic gas analysis provides crucial, independent ver- 
ification that stellar dynamics can be used to determine 
reliable SMBH masses, and effectively rules out patho- 
logical orbital structure for the central part of NGC 4649. 

Our results provide strong support not only for the 
stellar kinematical results but also the single-phase, hy- 

We note that, technically we also require the entropy profile not 
to behave pathologically at small radii. However, since hydrostatic 
equilbrium requires dS/dr > 0, the most discrepant behaviour we 
might expect to see in S is a precipitous drop. In such a case, 
this would tend to reduce a central peak in the model temperature 
profile, making the evidence for a supermassive black hole even 
stronger. 



drostatic equilibrium approximation we adopted. Firstly, 
the X-ray image shows a remarkably symmetrical mor- 
phology; we do not confirm the evidence of significant 
asymmetries previously reported with shallower data 
dRandaU et"aL||2004t [Shurkin et al.] |2008| ). Secondly, the 
hydrostatic mass model actually hts the data remark- 
ably well, despite the complicated temperature and den- 
sity profiles. Thirdly, we find good agreement between 
Mbh determined via both optical and X-ray techniques. 
As our derived mass profile was extended to larger radii, 
we found excellent agreement with stellar mass profile 
within ~3 kpc. At even larger radii, the mass determi- 
nation from the kinematic s of globular cluster s agreed 
very well with our results (Bridges et al. 2006). Taken 
together, these indicate that the ISM WGd649 must 
be very close to hydrostatic equilibrium at all scales, or 
there is a remarkable conspiracy in the data. 



In a recent provocative paper, Diehl fc Statlerj (2007) 
analysed the X-ray morphology ot a sample or early- 
type galaxies and concluded that, since the majority 
show clear evidence of morphological disturbances in 
their cores, probably all elliptical galaxies are out of 
hydrostatic equilbrium. Quite apart from the obvi- 
ously false syllogism that, since there exists disturbed 
(non-hydrostatic) ISM in many early-type galaxies, all 
early-type galaxies (including those without clear dis- 
turbances) are non- hydrostatic, NGC 4649 represents a 
clear counter-example to this argument. This is unsur- 
prising, since cosmological simulations of galaxy clus- 
ters have shown that the existence of asymmetries typi- 
cally does not translate to mass errors of more than 25% 
QTsai et al.||l994| ; [Evrard et aL||1996| ; |Nagai et al.]|2007| ), 
even m systems that are manifestly more disturbed than 
NGC 4649. Comparing optical and X-ray mass estimates 
for two systems (NGC 1399 and M87) that are much 



more distur b ed th an NGC 4649, a recent study by Chu-^ 
razov et al.| ( ^007 ) reaches a similar conclusion. Hence, 



even when low-level asymmetrical features are eventually 
detected with higher quality data, as they must, there 
is no evidence from either theoretical or observational 
studies that this implies a systematic error in the de- 
rived mass much larger than the other systematic errors 
we addressed in this paper. 

With Chandra we were unable to resolve the sphere of 
influence (^^60 pc) of the SMBH on the gas although, 
since the emissivity is proportional to the square of the 
gas density (which rises steeply towards the centre), we 
were able to probe the gas density within a factor ~2-3 
of it. Provided hydrostatic equilibrium holds, resolving 
the sphere of influence is not strictly necessary to place 
constraints on the central mass if one takes care to com- 
pute reliable average densities and temperatures within 
each radial bin. Ideally, however. X-ray imaging spec- 
troscopy at scales smaller than this scale is desirable so 
that the SMBH mass is not determined primarily by a 
single data-point. With such data one could also search 
for morphological disturbances at the smallest scales, 
which would indicate local deviations from hydrostatic 
equilbrium. Unfortunately, such imaging is beyond the 
capabilities of Chandra and will have to wait until the era 
of higher spatial resolution X-ray instruments. Nonethe- 
less, within the central ^200 pc, the dynamical time- 
scale is only 8 x 10^ yr, as compared to a cooling time 
~ 2 X 10^ yr so that, unless the gas is being constantly 



11 



stirred up, we would expect it to relax quickly into a 
state approximating hydrostatic equilibrium. Although 
ISM disturbances due to A GN activity appe ar common 
in early- type galaxies (e.g. Jones et al. 2002| ), these are 
typically seen on large scales, comparable to any radio 
lobes fwhich , although very weak in N GC 4649, are on 
~kpc scales: Stanger fc Warwick 1986| ). Given that the 
dynamical time-scale is shortest in the central parts of 
the galaxy, NGC 4649 would have to be in a very unusual 
state for AGN activity to stir up the gas within only the 
central ~100 pc but leave the X-ray morphology outside 
this range undisturbed. Although the current X-ray data 
do not allow us completely to rule this out, turbulence 
or bulk moti on would provide additional support against 
gravity (e.g. Nagai et al. 2007 ), lowering the central tem- 
perature for a given mass and entropy, and leading to 
an underestimate of Mbh- Since a deviation from this 
assumption would lead to poorer agreement with the op- 
tical data, our results indicate that the gas is very close 
to hydrostatic. 

In addition to being close to hydrostatic equilibrium, 
our results indicate that the gas must also be close to sin- 
gle phase. We tested this explicitly by adding an extra 
vapec spectral component in each annulus. We found 
that the fit did not improve significantly (Ax^=9 for 
30 fewer degrees of freedom) and, excepting the inner- 
most bin, the extra component was negligible in each 
ann ulus. supporting the single-phase assumption. Al- 



thou ;h we cannot rule out multi-phase gas within the 
innermost ^^200 pc, it is difficult to understand why gen- 
uinely multi-phase gas would be confined only to this re- 
gion. Fitting a single temperature model to data where 
multi-temperature models are actually required i n gen- 
eral s hould lead us to underestimate the mass (e.g. Allen 
998), so, if anything, the effect of multi-phase gas would 
e to make our measured Mbh more discrepant with the 
optical data (implying the gas is probably single-phase). 
It is worth noting that our models do indicate a steeply- 
rising temperature profile within the innermost bin and 
in such cases, two-temperature spectral models are often 
statistically better descri ptions of th e data than single- 
temperature models (e.g. Buotc 2000). However, one can 
still obtain reliable mass profiles by interpreting a single- 
temperature fit to the data provided one takes care to 
compare its parameters to a suitably averaged mod el, as 
we did in the present work (Gastaldello et al. 2007). 

One of the intriguing results from our analysis is an 
explanation of the central temperature spike seen in this 
system. As we show in Fig ^ if we omit the cen- 
tral black hole, our model temperature profile is much 



less ( :entrally-peakc d than the real data. T his co nfirms 



the predictions of Brighcnti fc Mathews ~( 1^55 ), who 
suggested that such features should be prevalent when 
Chandra temperature profiles of elliptical galaxies host- 
ing massive black holes are investigated. This detection, 
therefore, opens the possibility of using the same tech- 
nique to measure the central black hole mass in other, 
morphologically relaxed, early-type galaxies. Since sig- 
nificant X-ray halos are most often associated with the 
most massive galaxies, this will make such a method par- 
ticularly sensitive to the most massive black holes. How- 
ever, even for a black hole as massive as ~ 3 x IO^Mq, the 
sphere of influence is only ^^50 pc, requiring precise tem- 
perature constraints on a comparable scale. Within the 



limitations of current instrumentation, this effectively re- 
stricts the analysis to very X-ray bright galaxies within 
~20 Mpc, severely limiting the number of systems in 
which we might hope to see it. With the improved spa- 
tial resolution promised by proposed fu ture X-ray satel 



lites. such as the Generation-X mission ( Windhorst et al 
2006), though, such measurements should become rou- 
tine. 

In practice, the presence of a central temperature spike 
is not actually a sufficient diagnostic for a massive cen- 
tral black hole. As is clear from Fig ^, a mild central 
increase in the temperature can arise as a consequence 
of the centrally-peaked stellar mass distribution. In fact, 
if the stellar mass was significantly more cuspy (for ex- 
ample in a "powerlaw" galaxy, rather than the "cored" 
NGC 4649) , we would have seen a significant central tem- 
perature rise even without a black hole. More impor- 
tantly, the mass profile alone does not determine the 
shape of the temperature profile, but instead it is the in- 
terplay of the mass and entropy profiles. The shape of the 
latter is set by the balance between heating and cooling, 
which is not completely understood. In general, the flat- 
ter the central entropy profile, the more centrally-peaked 
is the temperature profile, so quantitative differences in 
the thermal histories of individual galaxies would explain 
in part the paucity of such temperature spikes. In many 
objects, this is compoun ded by the presence of AGN- 
driven distu rbances (e.g. Birzan et al. 2004; 'Sullivan 



et al. 



2005 ), which tend to stir up the gas, disturbing 



hydrostatic equilibrium and mixing hot, central gas with 
cooler gas from larger radii. Combined with the obser- 
vational difficulty of the measurement, these effects may 
expl ain the few reported galaxie s exhibiting this feature 
(e.g. Eumphrey et al. 2006 , 2004 ), in contrast to the ubiq- 
uity of black holes. 

In addition to the black hole mass estimate, our 
method also allowed us to determine directly the total 
stellar mass of the system, without making any assump- 
tions regarding the stellar mass-to-light ratio, and re- 
quiring only that stellar mass follows light. We therefore 
obtained both data-points (albeit not completely inde- 
pendently) to place on a Mbh versus stellar mass re- 
lation. We estimated that the mass of the SMBH is 
~1.7% times the mass in stars of the galaxy. It is in- 
teresting to compare this to the relations between Mr h 
and bulge mass (Mbui) found by [Marconi fc Hunt] ( |2003| ). 
For a bulge mass of ~ 2 x lO^^Mo, those authors would 
predict a ratio of Mbh to Mbui of ~0.2%, which is dis- 
crepant at almost an order of magnitude from our mea- 



suremen t (larger than the expected scatter). Marconi 



fc Hunt] estimated Mbui by inserting the stellar velocity 
dispersion and effective radius into the Virial relation, 
which actually gives an estimate of the total gravitating 
matter within, roughly, the effective radius. For the case 
of NGC 4649, which they include in their sample, they 
adopted an effective radius of 7.5 kpc (corrected to our 
adopted distance) and derived Mbui— 7.8 x lO^-^M©. This 
is considerably higher than the stellar mass we measured 
within the galaxy and, not unexpectedly, falls far short 
of the 3 x 10^^M(7, Virial mass of the entire system 



(Humphrey et al. 2006). Instead, this is approximately 
the total gravitating mass within ~19 kpc, of which dark 
matter comprises almost ~75%. Clearly this illustrates 
the difficulty in interpreting Virial theorem estimates for 
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Mbui, especially at the most massive end where galaxies 
sit in the centres of group-like dark matter halos. In or- 
der to gain insight into the SMBH-bulge relation, a full 
decomposition of the mass profile into stellar and dark 
components is, therefore, essential. 

Finally, based on the mean stellar population age 
a.nd abundance estimat ed from our Lick index analysis 
( [Humphrey et al.| |2006| ), we can estimate the expected 
mass-to-light (M/Lj) ratio for different simple stellar 
population models, assuming all stars were bo rn in a 



swi ft burst of star-formation. Using the models of [Maras- 



ton | \200^), the 

suming a Kroupa (2001) IMF, and depending only very 



J-band ratio is 1.86 ± O.21M0Lq , 



weakly on the mctallicity of the stars. As an alt erna- 



tive, if we adopt version 2 of the PEGASE code ( Fioc 
& R( )cca- Volmcrange 1997), the predicted J-band M/L 



ratio is 1.55±O.16M0Lq , indicating the level of system- 
atic uncertainty in such population models. These results 
compare favourably with our observed stellar M/Lj ra- 
tio (1.37 ± O.IOMqLq^), although the |Maraston| mod- 
els are slightly discrepant. We note that our measured 
M/L J is completely inconsistent with the predictions for 
a Salpeter IMF (2.7±0.3MoLq^ for both sets of models). 



as we found in Humphrey et al. ( 2006 ), unless the age 
for the stellar population is dramatically over-estimated. 
In that case it would need to be closer to ~7 Gyr than 
the 13±2 Gyr we measured. 
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